Nuclear suppression of dileptons at forward rapidities 
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O , Abstract. Data from E772 and E866 experiments on the Drell-Yan process exhibit a 



o 



significant nuclear suppression at large Feynman xf- We show that a corresponding kinematic 
region does not allow to interpret this as a manifestation of coherence or a Color Glass 
■ Condensate. We demonstrate, however, that this suppression can be treated alternatively as 

an effective energy loss proportional to initial energy. To eliminate suppression coming from 
the coherence, we perform predictions for nuclear effects also at large dilepton masses. Our 
calculations are in a good agreement with available data. Since the kinematic limit can be also 
approached in transverse momenta pr, we present in the RHIC energy range corresponding 
predictions for expected large-pT suppression as well. Since a new experiment E906 planned 
^ l' at FNAL will provide us with more precise data soon, we present also predictions for expected 

O , large-SF nuclear suppression in this kinematic region. 

1. Introduction 

^ ■ In comparison with a central region of very small rapidities, y — t- 0, the forward rapidity 

0> \ region allows to study processes corresponding to much higher initial energies accessible at 

' mid rapidities. If a particle with mass M and transverse momentum px is produced in a hard 

. reaction then the corresponding values of Bjorken variable in the beam and the target are 

^ ■ 3^1,2 = \J M"^ + pj, e^y/ y/s. Thus, at forward rapidities the target X2 is e^- times smaller than at 

mid rapidities. This allows to study already at RHIC coherence phenomena (shadowing. Color 
Glass Condensate (CGC)), which are expected to suppress particle yields. 

Forward rapidity physics, manifested itself as a strong nuclear suppression, has been already 
investigated in variety of processes at different energies: in production of different species of 
^ ■ particles in p + A collisions [1], in charge pion [2] and charmonium production [3l H] at SPS, 

in the Drell-Yan process and charmonium production at Fermilab [5l [6] and later on at larger 
RHIC energies by measurements of high-pj- particles in d -|- Au collisions [3 |8] . 

Althought forward rapidity region at RHIC allows to investigate small-x coherence 
phenomena, one should be carreful with interpretation of observed suppression. Such a 
suppression is arisen globally for any reaction studied so far at any energy. Namely, all fixed 
target experiments have too low energy for the onset of coherence effects since X2 is not small. 
The rise of suppression with y (with Feynman xp) shows the same pattern as observed at RHIC. 

This universality of suppression favors also another mechanism which should be common for 
all reactions studied at any energy. Such a mechanism was proposed in [9] and allows to describe 
a strong suppression via energy conservation effects in initial state parton rescatterings. It can 
be also interpreted alternatively as a parton effective energy loss proportional to initial energy 
leading so to scaling of nuclear effects. 

The projectile hadron can be decomposed over different Fock states. A nucleus has a higher 
resolution than a proton due to multiple interactions and so can resolve higher Fock components 
containing more constituents. Corresponding parton distributions fall off steeper at x — )• 1 



where any hard reaction can be treated as a large rapidity gap (LRG) process where no particle 
is produced within rapidity interval Ay = — ln(l — x). The suppression factor as a survival 
probability for LRG was estimated in [9], S{x) ~ 1 — x. Each of multiple interactions of 
projectile partons produces an extra S{x) and the weight factors are given by the AGK cutting 
rules [lOj. Then the effective parton distribution correlates with the nuclear target [9} lllj. 



fqfl{x,Q^,b) = C fg/N{x,Q'^)exp 



[l-S{x)]a,ffTA(b) 



(1) 



where Ta(6) is the nuclear thickness function defined at nuclear impact parameter b, CTeff = 
20 mb [U] and the normalization factor C is fixed by the Gottfried sum rule. 

In this paper we study a suppression of the Drell-Yan (DY) process on a nucleus with respect 
to a nucleon target and the rise of this suppression with y (xi, xp) in various kinematic regions. 
First we compare our predictions with data from the fixed target E772 experiment at FNAL [5] . 
Then similar nuclear effects are predicted also for the RHIC forward region expecting the same 
suppression pattern as seen at FNAL. Finally we perform for the first time predictions in the 
kinematic range corresponding to a new E906 experiment planned at FNAL where no coherence 
effects are expected. 

2. The color dipole approach 

The DY process in the target rest frame can be treated as radiation of a heavy photon/dilepton 
by a projectile quark. The transverse momentum pT distribution of photon bremsstrahlung in 
quark- nucleon interactions, a'^^ {a,'pT), reads [12]: 



da{qN ^^X) _ 1 r ,2 

in,f 



d{lna)d?pT (27r)2 



where S(a,ri,r2) = {(Tqq{ari) + aqq{ar2) - (Tqq{a\ri - f2|)}/2, a = pt^*lp\ and the light-cone 
(LC) wave functions of the projectile g+7* fluctuation <I>*.g(a, r) are presented in [12]. Feynman 
variable is given its, xp = x\ — X2 and in the target rest frame x\ = p'^* /pp. For the dipole cross 
section aqq{ar) in Eq. ([2]) we used GBW [13] and KST [T3] parametrizations. 

The hadron cross section is given convolving the parton cross section, Eq. ([2]), with the 
corresponding parton distribution functions (PDFs) fq and fq [12\ [15] , 



dxpd'^pxdM^ 3tt xi + X2 Jxi o"^ i a ' a' j d{lna)d'^p' 



Dp 



(3) 



where Zq is the fractional quark charge, PDFs fq and fq are used with the lowest order (LO) 
parametrization from [16j at the scale = + (1 — xi)M'^ and the factor aem/(3vrM^) 
accounts for decay of the photon into a dilepton. 

3. Dilepton production on nuclear targets 

The rest frame of the nucleus is very convenient for study of coherence effects. The dynamics 
of the DY process is controlled by the coherence length, 

^ _ 2Eq a(l - g) _ _i (1 - a) 

''~ {I- a)M'^ +a^ml+p\ ~ niN X2 {I - a) + ml + p^ ' ^' 

where Eq = Xqsj^mj^ and mq is the energy and mass of the projectile quark and the center of 
mass energy squared s = (M^ + Pp)/xiX2- The fraction of the proton momentum Xq carried by 
the quark is related to xi as axq = xi. 

The coherence length is related to the longitudinal momentum transfer, qp = 1/lc, which controls 
the interference between amplitudes of the hard reaction occurring on different nucleons. The 



condition for the onset of shadowing in a hard reaction is sufficiently long coherence length 
(LCL) in comparison with the nuclear radius, Ic > Ra, Here the special advantage of the color 
dipole approach allows to incorporate nuclear shadowing effects via a simple eikonalization of 
crqqix,r) yLZj, i-e. replacing cjqq{x,r) in Eq. ^ by cT^g(x,r): 

(5) 



'qq\ 
1 



1 -J^'^m^Aib) 



The corresponding predictions for nuclear broadening in DY reaction based on the theory [12] 
for LCL limit were presented in [18]. 

In the short coherence length (SCL) regime the coherence length is shorter than the mean 
internucleon spacing, Ic < 1-^2 fm. In this limit there is no shadowing due to very short duration 
of the 7* + fluctuation. The corresponding theory for description of the quark transverse 
momentum broadening can be found in [19', '20]. 

In this regime the transverse momentum distribution for an incident proton can be obtained 
integrating over a similarly as in Eq. Q: 

where a'^^{a,pT) represents the cross section for an incident quark to produce a photon on a 
nucleus A with transverse momentum pT- This cross section can be expressed convolving the 
probability function W'^^{kTiXq) with the cross section a'i^ {a^kx) (see Eq. ([2])), 

a^^(a,pT) = j d'kTW'i^{kT,Xq)a''^{a,rT) , (7) 

where It = Pt — ctk^- 

Probability distribution in Eq. ([7]) that a quark will acquire transverse momentum on the 
nucleus, W'^^{kT, Xq), is obtained by the averaging procedure over the nuclear density pA{b, z): 

W^^ikT. xq) = \j d^dzpAib, z) WlikT, Xq, b, z) , (8) 

where W\{kTi Xq, b, z) = duq/d'^kx means now the partial probability distribution that a valence 
quark arriving at the position (6, z) in the nucleus A will have acquired transverse momentum kj-- 
It can be written in term of the quark density matrix, ^q{ri,f2) = (&o/^) 6xp(— 6o(''i +1^2)/'^)^ 

Wl{kT,Xq,b,z) =777^72 / (iVidV2e*^^-(^^-"^2)^^(^^y2)e-^"«"'(^«'"^i-^2)^A(^+^'^), (9) 



where 6g = A , with (r^^) = 0.79 it 0.03 fm^ representing the mean-square charge radius of the 

^ ch' 

proton. Ta(6, z) in Eq. 1^ is the partial nuclear thickness function, Ta(6, z) = ^ dz' pA{b, z'). 

Nuclear effects in p+A collisions are usually investigated via the so called nuclear modification 
factor, defined as Ra{pt,xf,M) = f^J^al^^^M^ /A f^l^i^^dkfi , where the numerator is 
calculated in SCL and LCL regimes as described above. Corrections for the finite coherence 
length was realized by linear interpolation using nuclear longitudinal formfactor [2T] (for more 
sophisticate Green function method see [HI [22j ) . 

Note that at RHIC energy and at forward rapidities the eikonal formula for LCL regime, 
Eqs. ([3]) and ([5]), is not exact since higher Fock components containing gluons lead to additional 
corrections, called gluon shadowing (GS). The corresponding suppression factor Rq was derived 
in [211[18] and included in calculations replacing in Eq. ([5]) aqq by Rc^qq- GS leads to reduction 
of the Cronin effect [23] at medium-high pT and to additional suppression (see Fig. [3]) . 

In the fixed target FNAL energy range, for elimination of the coherence effects one can study 
production of dileptons at large M (see Eq. dH)) as has been realized by the E772 Collaboration 
|5]. Another possibility is to study the DY process at large xi — )• 1, when also a — )• 1, and /c — ^ 
in this limit (see Eq. (jl])). 
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Figure 1. Differential cross section of 
dileptons in p + p collisions at xp = 0.63 
and M = 4.8 GeV vs. E866 data 



4. Nuclear suppression at forward rapidities: model vs. data 

We start with the DY process in p + p collisions. Besides calculations based on Eq. ([3|) using 
GRV98 PDFs [16] (see the solid line in Fig. [1]) we present by the dashed and dotted line also 
predictions using proton structure functions from |25j and CTEQ6L parametrization of PDFs 
from [26], respectively. Fig. [1] shows a reasonable agreement of the model with data from the 
E866/NuSea Collaboration |24|. This encourages us to apply the color dipole approach to nuclear 
targets as well. 




Figure 2. Ratio R^^{W/D) of Drell-Yan cross sections on W and D vs. E772 data for 
6 < M < 7 GeV (Left) and 7 < M < 8 GeV (Right). Predictions correspond to SCL (dotted 
curves), LCL (dashed curves) regimes and their interpolation (solid curves). Thick and thin 
curves are calculated with and without corrections ([T]) for energy conservation, respectively. 



The E772 Collaboration [5j found a significant suppression of DY pairs at large xi (see 
Fig. [2]) . Large invariant masses of the photon allows to minimize shadowing effects (see a small 
differences between dotted and solid lines in Fig. [2]). If effects of energy conservation, Eq. ([1]), 
are not included one can not describe a strong suppression at large xi. In the opposite case a 
reasonable agreement of our model with data is achieved. 



One can approach the kinematic hmit increasing px- Therefore we present also predictions 
for pt dependence of the nuclear modification factor Rd+Au at RHIC energy and at several fixed 
values xp. Similarly as in [9] instead of usual Cronin enhancement, a suppression is found 
(see Fig. [3|). The onset of isotopic effects (IE) in d + Au collisions at large pr gives the values 
RdfAu ~ 0-73 ^ 0.79 depending on xp- In p + Au collisions the corresponding ratio Rp+Au — ^ 1 
from above and no nuclear effects are assumed at large pT expecting so QCD factorization. 
However, we predict a strong onset of effective energy loss effects, Eq. ([T]), at large xp (see 
Fig. [3D quantifying itself as a large deviation of suppression from the above values R^^j^w The 
predicted huge rise of suppression with xp in Fig. [3] reflects much smaller survival probability 
S{xp) at larger xp and can be tested in the future by the new data from RHIC. Note that effects 
of GS depicted in Fig. Elby the thick lines lead to additional suppression which rises with xp. 
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Figure 3. (Left) Predictions for the ratio Rd+AuiPr) at ^/s = 200 GeV for several fixed values 
xp without effects of effective energy loss. (Right) The same as (Left) but with effects of 
effective energy loss, Eq. (thin lines). Thick lines additionally include GS effects. 

Finally we present in Fig. [5] for the first time predictions for xi dependence of the nucleus- 
to-nucleon ratio in the kinematic range corresponding to a new E906 experiment planned at 
Fermilab. We shoud not expect any shadowing effects since initial energy is small, Eiab = 
120 GeV and a strong nuclear suppression at large xi is caused predominantly by the energy 
conservation constraints, Eq. ([1]). 

5. Summary 

We demonstrate that besides an onset of coherence a nuclear suppression at forward rapidities 
(large xi, xp) can be induced also by energy conservation effects in multiple parton rescatterings 
interpreted alternatively as a parton effective energy loss proportional to initial energy. 
Universality of this treatment is in its applicability to any reaction studied at any energy also 
in the kinematic regions where coherence phenomena (shadowing, CGC) can not be manifested. 
First we apply this approach to the DY process and explain well a significant suppression at large 
xi in accordance with the E772 data. The FNAL energy range and large invariant masses of 
the photon allow to minimize the effects of coherence, what does not leave much room for other 
mechanisms, such as CGC. Then we predict a significant suppression also for d + Au collisions 
at RHIC in the forward region (see Fig. [3]) . At small pT we show an importance of GS effects 
and their rise with xp. Finally we present for the first time predictions for strong nuclear effects 
expected in a new E906 experiment planned at FNAL. Much smaller beam energy than in E772 
experiment allows to exclude safely interpretations based on coherence phenomena. 
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Figure 4. (Left) Predictions for the ratio R^^{W/D) of Drell-Yan cross sections on W and D 
for 6 < M < 7 GeV (Left) and 7 < M < 8 GeV (Right) reahzed for the kinematic range of the 
planned E906 experiment at Fermilab. SoUd and dashed curves are calculated with and without 
effects of effective energy loss, Eq. ([T]), respectively. 
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